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Abstract. Star clusters are observed to form in a highly compact state and 
with low star-formation efficiencies, and only fO per cent of all clusters appear 
Oh to survive to middle- and old-dynamical age. If the residual gas is expelled on 

a dynamical time the clusters disrupt. Massive clusters may then feed a hot 
kinematical stellar component into their host-galaxy's field population thereby 
■ thickening galactic disks, a process that theories of galaxy formation and evolu- 

tion need to accommodate. If the gas-evacuation time-scale depends on cluster 
mass, then a power-law embedded-cluster mass function may transform within 
a few dozen Myr to a mass function with a turnover near 10 5 M , thereby possi- 
^ ' bly explaining this universal empirical feature. Discordant empirical evidence on 

the mass function of star clusters leads to the insight that the physical processes 
shaping early cluster evolution remain an issue of cutting-edge research. 
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1. Early cluster evolution 

The star- formation efficiency (sfe), e = M ec \/ (M cc \ + M gas ), where M ec \, M gas are 
the mass in freshly formed stars and residual gas, respectively, is 0.2 ^ e ^ 0.4 
£j, ! ( Lada Ladal 120031 ) implying that the physics dominating the star- formation 



process on scales < 10 pc is stellar feedback. Within this volume, the pre-cluster 
cloud core contracts under self gravity thereby forming stars ever more vigor- 
ously, until feedback ener gy suffices to halt the process (feedback-termination) 



( Weidner &: Kroupal20Qq ) . This occurs on one to a few crossing times (~ 10 6 yr) , 



and since e ach proto-star needs about 10 5 yr to accumulate about 95 per cent 



of its mass ( Wuchterl Tscharnuterll2003f ). the assumption may be made that 



the embedded cluster is mostly virialised at feedback-termination. Its stellar 
velocity dispersion, 



a^^GM ecl /(eR), (1) 

may then reach a = 40pc/Myr if M ec \ = 10 5 ' 5 Mq which is the case for eR < 
1 pc. This is easily achieved since the radius of one-Myr old clust ers is R ~ 1 pc 
with a weak, if any dependence on mass ((Bastian et al 1 l200^ . Very young 



clusters (age ^ 10 Myr) would thus appear super-virial, i.e. with a velocity 
dispersion too large for the cluster mass. 

The above exercise demonstrates that the possibility may be given that a 
hot kinematical component could add to a galactic disk as a result of clustered 
star formation for reasonable physical parameters. Thickened galactic disks may 
result. But this depends on 
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(ii) i? (cluster concentration) and 

(iii) the ratio of the gas-expulsion time-scale to the dynamical time of the 
embedded cluster, r gas /i cross . 

1.1. Empirical constraints 

The first (i) of these is clearly fulfilled: e < 40 per cent (lLada fc Ladal2003h . The 



second (ii) also appears to be fulfilled such that clusters with ages ^ 1 Myr have 
R ^ 1 pc ind ependen t ly of t heir mass. Some well-studied cases are tabulated and 
discussed in IKroupal (|2005h . Finally, the ratio r gas /t cross (iii) remains uncertain 
but critical. 

The well-observed cases discussed in IKroupal (|2005h do indicate that the 
removal of most of the residual gas does occur within a cluster-dynamical time, 
TgasAcross ~ 1- Examples noted are the Orion Nebula Cluster (ONC) and R136 
in the LMC both having significant super-virial velocity dispersions. Other ex- 
amples are the Treasure-Chest cluster and the very young star-bursting clusters 
in the massively-interacting Antennae galaxy which appear to have HII regions 
expanding at velocities such that the cluster volume may be evacuated within a 
cluster dynamical time. 

A simple calculation of the amount of energy deposited by an O star 
within a cluster crossing time into its surroundi ng cluster-nebu la also suggests 
it to be larger than the n ebula binding energy (|Kroupall2nnl '>. Furthermore, 
iBastian fc Goodwinl (l2006T > note that many young clusters have a radial-density 
profile signature as expected if they are expanding rapidly. 

Thus, the data suggest the ratio T gas /t cross to be near one, but much more 
observational work needs to be done to constrain this number. Measuring the 
kinematics in very young clusters would be an extremely important undertaking, 
because the implications of r gas /t cross ^ 1 are dramatic. 

To demonstrate these implications it is now assumed that a cluster is born 
in a very co mpact state (R « 1 pc), with a low sfe (e < 0.4) and T gas /£ cr0 ss ~ 1- 
As noted in IKroupal ((20051 1. "in the presence of O stars, explosive gas expulsion 
may drive early cluster evolution independently of cluster mass" . 



2. Implications 

2.1. Cluster evolution and the thickening of galactic disks 

As one of the important implications, a cluster in the age range of ~ 1 — 50 Myr 
will have an unphysical M/L ratio because i t is out of dynamical equ ilibrium 
rather than having an abnormal stellar IMF ( Bastian Goodwinl |2b06). 



Another implication would be that a Pleiades-like open cluster would have 
been born in a very dense ONC-type configuration and that, as it evolves, 
a "moving-group-F is established during the first few dozen Myr comprising 
roughly 2/3rds of the initial stellar population and expanding outwards with 
a velocity dispersion which is a f unction of the pre-gas-expulsion configuration 
( Krouoa. Aarseth &: Hurlevll20 01'). These computations were in fact the first to 



demonstrate, using high-precision A r -body modelling, that the re-distribution of 



3 



energy within the cluster during the embedded phase and the expansion phase 
leads to the formation of a substantial remnant cluster despite the inclusion of 
all physical effects that are disadvantageous for this to happen (explosive gas 
expulsion, Galactic tidal field and mass loss from stellar evolution). 

Thus, in this scenario stars form in very compact clusters that have radii 
less than about 1 pc and masses larger than a dozen Mq, and in the presence 
of O stars the residual gas is removed explosively leading to loss of stars from 
the cluster which form a moving group I (an expanding population of sibling 
stars). A cluster re-forms after gas expulsion as a result of energy-equipartition 
during the embedded phase, or during the expanding phase, and fills its Roche 
lobe in the Galactic tidal field. Such a cluster appears with an expanded core 
radius (compare the ONC with about 0.2 pc and the Pl eiades with 1.3 pc) and 
evolves secularly through evaporation until it's demise ( Baumgardt &: Makind 
2003; L amers et al.ll2005al lbh. A ' 'moving-group-II" establishes during this stage 
as the "classical" moving group made-up of stars which slowly diffuse/evaporate 
out of the re-virialised cluster remnant with relative kinetic energy close to 
zero. Under unfavorable conditions no cluster re-forms after gas expulsion. 
Such conditions may arise from time-variable tidal fields (e.g. through nearby 
cloud for mation) or simply through unfavorable density profiles just before gas- 
expulsion ((Boilv Kroupall2003al lbh . and they may dominate over the favorable 
ones (§ I2.3.|) . Note that the number of stars in moving-group-I always outnum- 
ber the number of stars in moving-group-II, while the unfavourable conditions 
would lead to no moving-group-II. 

Moving-groups-I would be populated by stars that carry the initial kinemat- 
ical state of the birth configuration into the field of a galaxy. Each generation 
of star clusters would, according to this picture, produce overlapping moving- 
groups-I (and II), and the overall velocity dispersion of the new field population 
can be estimated by adding in quadrature all expanding populations. This in- 
volves an integral over the embedded-cluster mass function, £eci(M ec \), which 
descri bes the distribution of the stellar mass content of clu sters when they are 
born (IKroupall200zl. l200a). It is know n to be a power-law (|Lada fc Ladal 12003: 
iHunter et al . 2003; IZhang &: FallllT999l ^ The integral can be calculated for a first 
estimate of the effect. The result is that for reasonable upper cluster mass limits 
in the integral, M ec i ^ 10 5 Mq, the observed age-velocity dispersion relation of 
Galactic field stars in the solar neighbourhood can be reproduced. 

This theory can thus explain the "energy deficit" , namely that the observed 
kinematical heating of field stars with age cannot, until now, be explained by 
the diffusion of orbits in the Galacti c disk as a res ult of scattering on molecu- 
lar clouds, spiral arms and the bar (|.IenkinsHl992T l. Because the age-velocity- 
dispersion relation for Galactic field stars increases with stellar age, this notion 
can also be used to map the star-formation history of the Milky- Way disk by re- 
sorting to the observed correlation between the star-formation rate in a galaxy 
and the maximum star-cluster m ass born in the population of young clusters 
( Weidner. Krouoa Larsenll2004h . 

A possible cosmolo gically-relevant imp lication of this "popping-cluster" mo- 
del emerges as follows (jKrouoal 1200^ .12005): A thin galactic disk which expe- 
riences a significant star burst such that star clusters with masses ranging up 
to 10 6 Mq can form in the disk, will puff-up and obtain a thick-disk made-up 
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of the fast-moving velocity wings of the expanding clusters after residual gas 
expulsion. A notable example wher e this proces s may have lead to thick-disk 
formation is the galaxy UGC 1281 ([Mould! I2005T ). which, by virtue of its LSB 
character, co uld not have obtained its thick d i sk as a result of a merger with an- 
other galaxy. lElmegreen Melov Elmegreenl (|2006f ) observe edge-on disk galax- 
ies with a chain-like morphology in the HST ultra-deep field finding many to have 
thick disks. These appear to be derived from kpc-sized clumps of stars. These 
clumps are probably massive star-cluster complexes which form in perturbed 
or intrinsically unstable gas-rich and young galactic disks. The individual clus- 
ters i n such complexes would "pop" as in the theory above ((Fellhauer &: Kroupal 
l2005h . leading naturally to thick disks as a result of vigorous early star formation. 

2.2. Structuring the initial cluster mass function 

Another potentially important implication from this theory of the evolution 
of young clusters is that if the gas-expulsion time-scale and/or the sfe varies 
with initial (embedded) cluster mass, then an initially featureless power-law 
mass function of embedded clusters will rapidly evolve to one with peaks, dips 
and turnovers at "initial" cluster masses that characterize changes in the broad 
physics involved, such as the gas-evacuation time-scale. 

Note that the "embedded" cluster mass is the birth mass in stellar con- 
tent, while the "initial" cluster mass is the mass of stars the virialised post-gas- 
expulsion star cluster would have had if it had been born with a SFE of 100 per 
cent (e = 1). 

"Initial" cluster masses are derived by dynamically evolving observed clus- 
ters backwards in time to t = by taking into account stellar evolution and 
a galactic tidal field but not the birth process ("classica l" evolution tracks, eg. 
Portegies Z wart et al. I l200ll : iBaumgardt fc Makmoll2003l). 

To further quantify this issue, iKrouoa fc Boilvi (12002) assumed that the 
function 

Mid = /st Med, (2) 
exists, where M ec \ is as above, M- 1C \ is the "classical initial cluster mass" and 

/st = / St (Med). (3) 

Thus, for example, for the Pleiades, M c \ ~ 1000 M Q at the present time (age: 
about 100 Myr), and a classical initial model would place the initial cluster 
mass at M- lc \ sa 1500 M & by using standard iV-body calculations to quantify the 
secular evaporation of stars from a n initially bound and virialised "classical" 
cluster (|Portegies Zwart et alJl200lft . If, however, the sfe was 33 per cent and 
the gas-expulsion time-scale was comparable or shorter than the cluster dynam- 
ical time, then the Pleiades would have been born in a compact configuration 
resembling the QNC and with a mass of embedded stars of M ec \ w 4000 M & 
( Krouoa. Aarseth fc Hurtevll200ll l. Thus, / st (4000M o ) = 0.38 in this particular 
case. 

By postulating that there exist three basic types of embedded clusters, 
namely 
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• clusters without O stars (type I: M ec \ ^ 10 2,5 M , e.g. Taurus- Auriga pre- 
main sequence stellar groups, p Oph), 

• clusters with a few O stars (type II: 10 2 - 5 ^ M ecl /M© < 10 55 , e.g. the 
ONC), and 

• clusters with many O stars and with a velocity dispersion comparable to 
the sound velocity of ionized gas (type III: M ec \ ^ 10 5 ' 5 Mq), 

it can be argued that / st ~ 0.5 for type I, / s t < 0.5 for type II and / st ~ 0.5 for 
type III. 

The reason for the high / st values for types I and III is that gas expulsion 
from these clusters may be longer than the cluster dynamical time because 
there is no sufficient ionizing radiation for type I clusters, or the potential well 
is too deep for the ionized gas to leave (type III clusters). Type I clusters are 
excavated through the cumulative effects of stellar radiation and outflows, while 
type II clusters pro bably require multiple supernovae of type II to unbind the 
gas ( Goodwinlll997f ). According to the present notion, type II cl usters undergo 



a dis ruptive evolution and witness a high "infant mortality rate" (|Lada Sz Ladal 
l2003h . therewith being the pre-cursors of OB associations and open Galactic 
clusters. 

Under these conditions and an assumed inverted Gaussian functional form 
for / s t = fst(M ec \), the power-law embedded cluster mass function transforms 
into a cl uster mass function wit h a turnover near 10 5 M and a sharp peak near 
1O 3 M (iKrouoa fc Boirvll200ifl . This for m is strongly r eminiscent of the initia l 
globular cluster mass functio n which is inferred by e.g. IVesperinl (|l998l . [2001); 
iParmentier &: Gilmord (2005); Baumgardt & Makino (2003) to be required for 



a match with the evolved cluster mass function that is seen to have a universal 
turnover near 10 5 Mq. Thus, ~ 10 10 yr-old cluster populations have log- normal 
Gaussian mass functions quite independently of the environment. 

In the theory presented here, the power-law embedded cluster mass func- 
tion evolves in the first 10 7-8 yr to the "initial" cluster mass function with 
the broad turn-over near 1O 5 M . In this context it is of interest to note that 
Ide Griis et aD (l2003h found a turnover in the cluster mass function near 10 5 Mq 
for clusters with an age near 10 9 yr in the star-burst galaxy M82. This result 
was challenged, but in their section 2.1 Ide Griis et all ((20 05) counter the criti- 
cisms. This age of 1 Gyr is too young to be explainable through secular evolution 
starting from a power-law cluster mass function. In contrast to this interpreta- 
tion, the theory outlined here would allow a power-law embedded cluster mass 
function to evolve to the observed form within 1 Gyr. Evaporation through two- 
body relaxation in a time-variable tidal field then takes over th ereby enhancing 
the observed log- normal shape of t he old cluster mass function (IVesDerinil H998. 
120011 : IParmentier k Gilmoreil200l eg.). 

This ansatz may thus bear the solution to the long-standing problem that 
the initial cluster mass function needs to have this turnover, while the observed 
mass functions of very young clusters are featureless power-law distributions. 

2.3. Cluster disruption independently of mass 

If, on the other hand, / s t = constant such that only about 10 per cent of all clus- 
ters survive independently of mass, then the above long-standing problem would 
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remain unsolved. This situation is suggested by the results of lLada fc Ladal 
( 2003) and Fall et al. I (|2005l ). among others. Observations show that 90 per cent 
of all very young clusters dissolve within about 10 yr, independent ly of their 
mass. This may either be due to a constant / st , or, as pointed out bv lFall et al.l 
( 20051 ). a result of 90 per cent of all clusters dissolving and only the rest surviving 
intact. The former case would simply imply a constant shift of the power-law 
mass function along the mass axis to smaller masses as a result of cluster dis- 
ruption through stellar feedback (operating in all embedded clusters), while the 
latter would imply a purely stochastic element to the disruption of clusters. 

It would therefore appear that this 'high infant mortality of clusters inde- 
pendently of c l uster mass 1 scenario would be in contradiction with the results of 
Ide Griis et alJ (2005) on the M82 1 Gyr-old cluster system combined with the 
observed power-law mass function of very young clusters. 



3. Conclusions 

Observations show that star clusters of any mass are formed extremely com- 
pact ( ^ 1 pc), while older clusters appear well distended with radii of a few-to- 
many pc. Very young clusters seem to be super-virial, and expanding HII shells 
also indicate the explosive removal of residual gas. Star-cluster formation thus 
appears to be rather violent leading to a loss of a large fraction of stars and 
probably in many cases total cluster disruption. 

As this discussion shows, the formation process of star clusters may have 
significant cosmological implications in that the morphology of galaxies may 
be shaped by the clusters born within them: the disks of galaxies are proba- 
bly thickened during periods of high star-formation rates while low-mass dwarf 
galaxies may attain halos of stars as a result of the violent processes associated 
with cluster birth as a result of which stars may spread outwards with relatively 
high velocity dispersions. Particular examples of such events have been noted 
above. If this is true can be verified by measuring the velocity dispersions of 
stars in clusters younger than at most a few Myr. This is very important to do 
because thickened galactic disks are often taken to be evidence for the merging 
of cold-dark-matter sub-structures. Naturally, in-falling dwarf galaxies will lead 
to thickened disks, and counter-rotating thick disks would be prime examples 
of such processes (|Yoachim & Dalcantonl (2005), but this discourse shows that 



sub-pc-scale processes would probably need to be incorporated in any galaxy- 
evolution study. 

The mass function of old globular clusters has a near-universal turnover near 
10 5 Mq which has so far defied an explanation, given that the mass function of 
very young clusters is observed to be a power-law, and that the modelling of clas- 
sical disruption through two-body relaxation and tidal fields poses a challenge 
for evolving a power-law to the approximate log-normal form for old clusters. 
Again, mass-dependent processes related to the physics of residual gas expul- 
sion can explain such a turnover, but in this case the turnover would have to 
b ecome evident in clust er populations younger than about 100 My r. The results 
of Ide Griis et al.l |2005) suggest this to be the case. The results of lZhang fc Fall 
( 19991 ). lowever, appear to show no such evidence for the cluster population in 



the Antennae galaxies: the mass function for two age groups of clusters (younger 
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than about 10 Myr, and between 10 and 160 Myr) show the same power-law 
form, albeit with different normalisations as a result of the high infant mortal- 
ity rate. If this is true, the n orbital aniso t ropics in the young globular cluster 
populations as proposed by iFall fc Zhansl (|200lT l would be needed to evolve a 



power-law mass function to the log-normal form. 

It therefore appears that the broad principles of early cluster evolution 
remain q uiet unclear. Given that clusters are the fundamental building blocks of 
qalaxies (lKrounall2005h this is no satisfying state of affairs. On the positive side, 



the importance of star clusters to stellar populations and other cosmological 
issues (merging histories of galaxies, galactic morphology) means that much 
more theoretical and observational effort is needed to clarify which processes 
act during the early life of a star cluster, and how these affect its appearance at 
different times. 
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